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ABSTRACT 



We present observations of SDSSp J104433. 04-012502. 2, a luminous quasar 
at z = 5.80 discovered from Sloan Digital Sky Survey (SDSS) multicolor imaging 
data. This object was selected as an z'-band dropout object, with i* = 21.8 ±0.2, 
z* = 19.2 ±0.1. It has an absolute magnitude M 1450 = —27.2 (H = 50 km s _1 
Mpc -1 , go — 0.5). The spectrum shows a strong and broad Lya emission line, 
strong Lya forest absorption lines with a mean continuum decrement Da = 0.91, 
and a Lyman Limit System at z — 5.72. The spectrum also shows strong 01 
and SilV emission lines similar to those of quasars at z ^ 5, suggesting that 
these metals were produced at redshift beyond six. The lack of a Gunn-Peterson 
trough in the spectrum indicates that the universe is already highly ionized at 
z ~ 5.8. Using a high-resolution spectrum in the Lya forest region, we place a 
conservative upper limit of the optical depth due to the Gunn-Peterson effect of 
r < 0.5 in regions of minimum absorption. The Lya forest absorption in this 
object is much stronger than that in quasars at z ^ 5. The object is unresolved 
in a deep image with excellent seeing, implying that it is unlensed. The black 
hole mass of this quasar is ~ 3 x 1O 9 M if we assume that it is radiating at the 
Eddington luminosity and no lensing amplification, implying that it resides in 
a very massive dark matter halo. The discovery of one quasar at M1450 < —27 
in a survey area of 600 deg 2 is consistent with an extrapolation of the observed 
luminosity function at lower redshift. The abundance and evolution of such 
quasars can provide sensitive tests of models of quasar and galaxy formation. 



Introduction 



At what epoch did the first generation of galaxies and quasars form? How was the 
universe re- ionized, ending the "dark ages" ( |Rees 1998| )? These fundamental questions can 
only be answered with studies of high-redshift objects. The last few years have witnessed 



the first direct observations of galaxies at redshift higher than five (Pey et al. 1998 



Weymann et al. 1998|, |Spinrad et al. 19"98|, |Ohen, Lanzetta fc Pascarelle 1999|, |van Breugel 



et al. 199~£ , Hu et al. 1999 , see also the review by |Stern fc Spinrad 1999 ), while detailed 
studies of the ensemble properties and large scale distribution of galaxies at z ~ 4 have 
begun (Steidel et al. 1998, 1999). Several quasars have been found at z ^ 5 (Fan et al. 1999, 
2000a, |Zheng et al. 2000]) , including a low-luminosity quasar at z — 5.50 flStern et al. 2000] ). 
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Studies of high-redshift quasars provide important probes of this critical epoch in 
cosmic evolution. The lack of the Gunn-Peterson (1965) effect in the absorption spectrum of 
a z ~ 5.0 quasar ( [Songaila et al. 19"9~9] ) indicates that the universe is already highly ionized 
at that redshift. The exact epoch of re-ionization could be determined from the absorption 
spectra of quasars at even higher redshift flMiralda-Escude 19971 , [Haiman fc Loeb 1999|) . 



The study of the luminosity function of high-redshift quasars will constrain models of quasar 
and galaxy evolution QHaiman fc Loeb 1998| , |Haehnelt et al. 1999|) , and determine whether 



it was UV radiation from AGNs or from young massive stars that re-ionized the universe, 
ending the "dark ages" ( |Haiman fc Loeb 1998|) . Measurements of the chemical abundance 



in the quasar environment will reveal the metal production process at the very early stage 
of galaxy evolution ( [Hammann fc Ferland 1999]) . Finally, luminous high-redshift quasars 
represent high peaks of density fields, and may be the markers of large scale structure at 
these early epochs ( |Djorgovski et al. 1999| , [Haiman fc Hui 2000| , |Martini fc Weinberg 2000] ) . 



The Sloan Digital Sky Survey (SDSS; [York et al. 2000| ) is using a dedicated 2.5m 
telescope and a large format CCD camera (|Gunn et al. 1998]) at the Apache Point 



Observatory in New Mexico to obtain images in five broad bands (V, g', r', i' and z', 
centered at 3540, 4770, 6230, 7630 and 9130 A, respectively; [Fukugita et al. 199Q ) over 
10,000 deg 2 of high Galactic latitude sky. The multicolor data from SDSS have proven to be 
very effective in selecting high-redshift quasars: more than 50 quasars at z > 3.5 have been 
discovered to date from about 600 deg 2 of imaging data (Fan et al. 1999, 2000a, [Schneider 



et al. 2000] , |Zheng et al. 2000]) . The inclusion of the reddest band, z', in principle enables 



the detection of quasars up to z ~ 6.5 in SDSS data. 

In this paper, we report the discovery of SDSSp J104433. 04-012502. 2 (the name 
reflecting its J2000 coordinates from the preliminary SDSS astrometry, accurate to ~ 0.1" 
in each coordinate), a very luminous, "z'-dropout" quasar at z — 5.80, selected by its very 
red i* — z* color. In a A-dominated flat universe (Hq = 65 km s" 1 Mpc -1 , A = 0.65 and 
Q = 0.35, referred to as the A-model in this paper, |Ustriker fc Steinhardt 1995| , [Krauss <5 



Turner 1995| ), z = 5.80 corresponds to an age of 0.9 Gyr in an universe 13.9 Gyr old, or a 



look-back time of 93.2% of the age of the universe. Similarly, the universe was 0.7 Gyr old 
at z = 5.80 in a universe 13.0 Gyr old at present for a model with Q — 1 and Hq = 50 km 
s -1 Mpc -1 , which we refer to as the Q = 1 model in this paper. We present the photometric 
observations and target selection in §2, and the spectroscopic observations in §3. In §4, 
we discuss the cosmological implications, including the constraints on the Gunn-Peterson 
effect, quasar evolution models, and black hole formation. 
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2. Photometric Observation and Target Selection 

The object SDSSp J104433.04-012502.2 (hereafter SDSS 1044-0125 for brevity) was 
selected from the SDSS imaging data based on its extremely red i* — z* color. The 
photometry of this object is summarized in Table 1. The photometric observations of this 
region were obtained by the SDSS imaging camera on 2000 March 4 during the SDSS 
commissioning phase. For off-equator scans, the telescope moves along a great circle, with 
the photometric camera drift- scanning at the sidereal rate. The effective exposure time 
is 54.1 seconds in each band. The seeing in the i' and z' bands was about 1.9". The 
photometric calibration is provided by an auxiliary 20-inch telescope at the same site 
(Uomoto et al. 2000, in preparation). The photometric zeropoint is accurate to about 7% 
in u' and z' and 2% in g', r' and i'. Because the definition of the photometric system 
is not yet finalized, we quote measured magnitudes using asterisks (e.g., i*) to represent 
preliminary photometry, while referring to the filters with primes (e.g., i'). SDSS 1044-0125 
is undetected in u', g' and r'. The object is detected at the 5-a level in the %' band. The 
z' detection is of very high significance, with i* — z* = 2.58 ± 0.20. Data are quoted as 
asinh magnitudes ( [Lupton, Gunn fc Szalay 1999|) and are on the AB magnitude system 



( |Fukugita et al. 1996|) . Finding charts for SDSS 1044-0125 in the i' and z' bands are shown 
in Figure 1. 



As shown by |Fan (199"9~)| , the colors of high-redshift quasars are strong functions of 



redshift in the SDSS filter system, as first the Lya forest and then the Lyman Limit 
Systems move through the filter system. At z > 3.6, quasars become very red in g* — r* 
while remaining blue in r* — i*, and can be readily distinguished from stars based on these 
colors. At z > 4.6, quasars become very red in r* — i* . They are easily distinguished from 
red stars in the r* — i* vs. i* — z* color-color diagram (Fan et al. 1999, 2000a). Quasars 
at z ~ 5.5 have very similar r* — i* and i* — z* colors to those of late M stars. Finally, at 
z ^ 5.7, the Lya emission line begins to move out of the SDSS %' filter. With a predicted 
j* _ z * 2 ; quasars at such redshifts become i'-band dropout objects. However, unlike the 
lower-redshift g' and r'-dropout quasars, these quasars have only one measurable optical 
color, as they will be completely undetected in r' . The SDSS alone cannot provide a 
constraint on the quasar's continuum shape redward of the Lya emission. It is thus difficult 
to distinguish them from other classes of red objects without additional information, such 
as near-infrared photometry (see also Zheng et al. 2000) or detection in the radio or X-ray. 

In fact, the only other known class of stellar objects with such red colors {i* — z* ^ 2) 
are the extremely cool stars/substellar objects with spectral type L or T ( |Kirkpatrick et] 



al. 1999| , |Strauss et al. Tggg , ll'svetanov et al. 2000| |Fan et al. 2000T| , |Leggett et al. 2U0q ) 



Most of these objects are brown dwarfs with mass below the hydrogen burning limit. 
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Although L and T dwarfs are very rare on the sky (the density for L dwarfs is 1 per 
~ 15 deg 2 for i* < 20, [Fan et al. 2000 b]) , they are still more numerous than are z > 5 



quasars, and are the major contaminants of searches for z'-dropout quasars. Near infrared 
photometry can be used to separate high-redshift quasar candidates from these cool dwarfs: 
the continuum shape of quasars is relatively flat towards near-IR bands, while the flux of L 
dwarfs continues to rise sharply towards longer wavelengths. 

We have selected i'-dropout candidates with i* — z* ^ 2 and z* < 19.5 from about 
600 deg 2 of SDSS imaging data. These 600 deg 2 overlap the publicly released area of the 
Two Micron All Sky Survey (2MASS). SDSS 1044-0125 is the only i'-dropout source in the 
2MASS covered area that is not detected at the 7-a level in any band in the 2MASS Point 
Source Catalog. All other sources are detected in at least one band at more than 10 a (e.g. 
''an et al. 20(Ub], |Leggett et al. 2000|). 



Is the non-detection in 2MASS sufficient to rule out the possibility that SDSS 
1044-0125 is a brown dwarf? In this area of the sky, the 2MASS Point Source Catalog has 
7-a limiting magnitudes of roughly J = 16.7, H = 15.9 and K s = 15.0. This indicates a 
7-a upper limit on the optical-IR colors of SDSS 1044-0125: z* - J < 2.5, z* - H < 3.3 
and z* — K s < 4.2. There are three L dwarfs in Fan et al. (2000b) that have i* — z* > 2: 
SDSS 0330-0025 (i* - z* = 2.13, spectral type L2), SDSS 0539-0059 (z* - z* = 2.31, L5) 
and SDSS 1326-0038 (i* - z* = 2.61, L8). They all have z* - J £ 2.7, z* - H £ 3.6 and 
z* — K s ^ 4.2 (see Tables 2 and 5 in [Fan et al. 2000b| ). Although the current L dwarf sample 



is small, the i* — z* and the z* — J (or H, K s ) colors seem to be reasonable indicators of 
the spectral type for L dwarfs ([Kirkpatrick et al. 1999[ |Fan et al. 2000b|) Therefore, if SDSS 



1044-0125 (i* — z* = 2.6) had been an L dwarf, we would expect it to have been detected by 
2MASS with high significance. Thus the non-detection in the 2MASS passbands indicates 
that SDSS 1044-0125 is unlikely to be a L dwarf, but is rather a source with a much flatter 
IR continuum, such as a quasar at z > 5.6 or a compact galaxy at z > 1. 

A K'-band image of SDSS 1044-0125 was obtained on the night of 2000 April 17, 
using the Near Infrared Camera (NIRC, Matthews & Soifer 1994) on the Keck I telescope 
under photometric skies and good seeing. The observations consisted of a nine-point 
dither pattern, integrating for 3 x 10s coadds at each location. The data were flattened, 
sky-subtracted, shifted, and stacked using the DIMSUM package in IRAF. The images 
show SDSS 1044-0125 to be an unresolved point source with FWHM = 0'.'375. Photometry 
through a 2" radius aperture yields K' = 17.02 ± 0.04, referenced to the standards of 
Persson et al. (1998). There are no companions or associated structure within 20" of the 
quasar to a 3a point source limiting magnitude of K' ~ 22.3. 
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3. Spectroscopy 

Low and high dispersion spectra of SDSS 1044-0125 were obtained with the Echelle 
Spectrograph and Imager (ESI; |Epps fc Miller 1998j ) on the Keck II telescope on the night 
of 2000 April 6. The night was photometric with 0.9" seeing. A 1200 second exposure was 
taken through a 0.7" slit in the low-dispersion mode of ESI. In this mode, a prism is used for 
dispersion. The spectrum covers 3900 A to about 10000 A. The dispersion varies roughly 
linearly with wavelength from 0.8 A/pixel at 3900 A to 10 A/pixel at 10000 A. In addition, 
two 1200 second high resolution spectra were taken through a 1.0" slit in the echellette 
mode of ESI. In this mode, the spectral range of 3900 A to 11000 A is covered in ten 
spectral orders with a constant dispersion of 11.4 km s _1 pixel -1 . Wavelength calibrations 
were performed with observations of Hg-Ne-Xe lamps in the low-dispersion prism mode, and 
a Cu-Ar lamp in the echellette mode. The spectrophotometric standard G191-B2B ( [Massey 



1988, Massey ik, Gronwall 1990| ) was observed for flux calibration. All observations were 



carried out at the parallactic angle. The data were reduced with standard IRAF routines. 

Figure 2 shows the final spectrum combining the low and high dispersion observations 
(total exposure time of 3600s), over a wavelength range of 4500 - 10000 A, and binned to 
4 A/pixel. The telluric absorption bands were removed using the standard star spectrum. 
The absolute flux scale of the spectrum is adjusted so that it reproduces the SDSS z* 
magnitude. The signal-to-noise ratio at A > 8000 A is 15 - 20 per pixel. The spectrum 
of SDSS 1044-0125 shows the unambiguous signature of a very high redshift quasar: the 
broad and strong Lya+NV emission line at A ~ 8300 A with a sharp discontinuity to the 
blue side, due to the onset of very strong Lya forest absorption. The flux level drops by a 
factor of ~ 10 from the red side to the blue side of the Lya+NV emission. A Ly/3+OVI 
emission line is detected at ~ 7000 A, with an additional flux decrement due to the onset 
of Ly/3 forest absorption lines. The spectrum shows no detectable flux at A < 6100 A 
because of the presence of a Lyman Limit System (see §4). Redward of Lya, two additional 
emission lines, OI+SiIIA1302 and SiIV+OIV]A1400 are also clearly visible. The synthetic 
i* — z* color calculated from the spectrum in Figure 2 is 2.5, consistent with the SDSS 
measurement. 

The redshift determination of such an object is not straightforward. The Lya emission 
line is severely affected by the strong Lya forest lines (as well as by possible internal 
absorption lines). For z ~ 4 quasars, Schneider, Schmidt & Gunn (1991) show that the 
peak of the Lya emission line is typically at rest-frame 1219 A, but with a large scatter. 
Another strong line, CIVA 1549, is out of the range of our spectrum. OI+SiIIA1302 and 
SiIV+OIV]A1400 are relatively weak. We use the central wavelengths of these two blends 
from [Francis et al. ( 1991 )| . Gaussian fits of these two lines are listed in Table 2. Based on 
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these fits, we adopt a redshift of 5.80 ± 0.02 for SDSS 1044-0125. Objective measurement 
of the equivalent width of the Lya line is difficult (see the discussion by Schneider, Schmidt 
& Gunn 1991). Using the continuum level determined from the red side of Lya emission, 
we find the rest-frame equivalent width of Lya +NV is ~ 26 A. However, this number is 
clearly an underestimate due to the extreme Lya absorption, maybe by a factor of two. 
Figure 2 indicates that the expected line center of Lya at z — 5.80 is severely absorbed, 
and that the blue wing of the line has almost completely disappeared. The measurement of 
the Ly/3+OVI line is difficult as well. We smooth the spectrum of Figure 2 to estimate the 
continuum level, and find an equivalent width of ~ 30A. This value is also highly uncertain 
due to the difficulty in defining a "continuum" in the Lya forest region. 

It is interesting to note that the strengths of the OI+SiIIA1302 and SiIV+OIV]A1400 
line blends are comparable to those at much lower redshift: the average rest frame equivalent 
widths of these two lines are 3.2 ± 0.4 A and 8.1 ± 0.6 A in the sample of 30 quasars at 
z ~ 4 in Schneider, Schmidt & Gunn (1991), compared to 3.4 and 7.0 A in the case of SDSS 
1044-0125. Previous studies have shown that quasar environments at z ~ 4 have roughly 
solar or higher metallicities ([Mammann fc Ferland 1999| ). Although we cannot derive the 
metallicity of SDSS 1044-0125 based on these emission line measurements, the existence 
of strong lines suggests that in this system, the metallicity is already quite high at this 
redshift. Assuming the metals were produced in stellar evolution, the initial starburst and 
chemical enrichment of the quasar environment must have happened at a very early epoch. 

From the spectrum in Figure 2, we derive its continuum AB magnitude at rest frame 
1280 A, ABi2m = 19.28, after correcting for interstellar extinction (E(B - V) = 0.054, 
from the map of [Schlegel, Finkbeiner fc Davis 199~8| ). In the A-model (see §1), SDSS 



1044-0125 has an absolute magnitude M 12 so = —27.41. Assuming a power law continuum 
with f v oc z/-°- 5 we find M 1450 = -27.50 and M B = -27.96. In the Q = 1 model (§1) it 
has M1280 = —27.15, Mi 45 o = —27.24, and Mb = —27.70. In this cosmology, the nearby 
luminous quasar 3C273 has Mb = —27.0. SDSS 1044-0125 is a very luminous quasar, about 
twice as luminous as 3C 273 (assuming that it is not amplified by lensing or beaming). 



SDSS 1044-0125 is detected neither in the FIRST radio survey ( [Becker, White fc 



Helfand 1995| ) at the lmJy level at 20cm wavelength, nor in the ROSAT All Sky Survey 



( |Voges et al. 1999]) , implying a 3-a upper limit of 3 x 10 13 ergs cm 2 s 1 in the 0.1 - 2.4 



keV band. These result is not unexpected; only a few z > 4 quasars have observed X-ray 
or radio fluxes above this value ( [Kaspi, Brandt fc Schneider 2000| , [Schmidt et al. 199H| ). A 



deep exposure with Chandra or XMM is needed to determine its X-ray properties. 
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4. Discussion 

4.1. Absorption Properties and Gunn-Peterson Effect 

The high luminosity of SDSS 1044-0125 makes it an ideal object for high signal-to-noise 
ratio observations to study the intergalactic medium at high redshift. In order to detect 
continuum break caused by the Lyman Limit System, an edge filter with width of 40 A 
was convolved with the low resolution spectrum in Figure 2. A strong peak at 6131 A is 
detected in the convolved spectrum, indicating the existence of a Lyman Limit System at 
zlls = 5.72. No flux is detected blueward of this break. A Lyman Limit System is usually 
detected within 0.1 of the emission redshift in essentially all quasars at z > 4 ( |Schneider^| 
Schmidt fc Gunn 1991| , |Storrie-Lombardi et al. 1996| , |Fan et al. 1999j ). 



The most striking feature of the spectrum of SDSS 1044-0125 is the very strong 
absorption caused by Lya forest lines. However, the flux level in the Lya forest region never 
reaches zero. It lacks the Gunn-Peterson (1965) trough that would exist in the spectrum 
of a quasar at redshift higher than the re- ionization redshift ([Haiman fc Loeb 1999 ), 



indicating that the intergalactic medium is already highly ionized at z ~ 5.8. We estimate 
the average continuum decrements as: Da.b = (1 — f° bs / f™/i where f° bs and fl° n are the 
observed and the unabsorbed continuum fluxes of the quasar, and Da and Db measure the 
decrements in the region between rest-frame Lya and Ly/3 (A = 1050 — 1170 A) and between 
Ly/3 and the Lyman Limit (A = 920 - 1050 A), respectively QOke fc Korycansky 19821) . The 
measurements of Da and Db require knowledge of the continuum shape redward of Lya. 
However, with Da approaching unity, the effect of different slopes is quite small. Assuming 
a power law continuum v a with a = —0.5, as indicated in Figure 3, we obtain Da = 0.91 
and Db = 0.94. Using a slope of -1.0 only changes the Da and Db values to 0.92 and 0.96, 
respectively. We therefore adopt Da = 0.91 ± 0.02 and Db = 0.95 ± 0.02. These values 
are in close accordance with those from quasar RD J030117+002025 (z = 5.50, |Stern et 



al. 2000| ) and from distant galaxies in the Hubble Deep Field ( |Weymann et al. 1998| ), and 



are much higher than that of the z = 5.00 quasar SDSSp J033829.31+002156.3 (D A = 0.75, 
Songaila et al. 1999| ), suggesting that the strong evolution of the strength of the Lya forest 



at z > 5, N(z) oc (1 + z) 2 ' 3 2 75 , measured at redshifts below five continues to a redshift of 
nearly six. Assuming this number density evolution, [Zuo (1993)1 an d |Fan (1999]] show that 



at z ~ 5.8, the expected average Da ranges from 0.8 to 0.9. 

We further derive an upper limit on the Gunn-Peterson optical depth following Songaila 
et al. (1999). Figure 3 shows the high-resolution echellette spectrum of SDSS 1044-0125 
in the Lya forest region (binned to 2 A/pixel). The continuum level is approximated by a 
jy-o.5 p OW6r i aw as a bove. Even the most transparent part of the forest does not return close 
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to the continuum level at this resolution. It is evident that the Lya forest is much stronger 
in SDSS 1044-0125 than in SDSSp J033829.31+002156.3 (Figure 2 of |Songaila et al. 1999| ), 
where a fraction of the forest has flux comparable to the extrapolated continuum. In the 
region between 7926 A and 7929 A, SDSS 1044-0125 has an optical depth r = 0.35 ± 0.07, 
where the error bar only reflects the statistical noise in the spectrum (note the noise 
level indicated in the figure). This value changes to 0.40 and 0.31 for power law slopes 
of 0.0 and —1.0, respectively. Therefore, we adopt a conservative limit of r < 0.5 at this 
redshift of 5.52. With higher resolution and signal-to-noise ratio, we might be able to select 
regions even less affected by the Lya forest lines. Therefore, it is only an upper limit. For 
comparison, Songaila et al. derived r < 0.1 for z = 4.72 with similar resolution. 

The Gunn-Peterson effect analysis above is based on an attempt to measure the 
amount of flux between Lya forest lines (e.g., |Giallongo et al. 1994j). At redshift higher 



than five, even with a moderately high resolution spectrum, these forest lines overlap, 
making it impossible to find a truly "line-free" region. Modern hydrodynamic simulations 
and semi-analytic models show that under the influence of gravity, the intergalactic medium 
becomes clumpy, and the Gunn-Peterson optical depth should vary even in the lowest 
column density regions (e.g. |Bi, Borner fc Ghu 1992| , Miralda-Escude fc Rees 1993| , |Gen| 
|et al. 1993| , [Hernquist et al. 1996| ). The minimum absorption regions in the forest merely 



represent regions that are most underdense in this fluctuating Gunn-Peterson effect. An 
accurate measurement of the Gunn-Peterson effect and the ionizing background from the 
high resolution spectrum of SDSS 1044-0125 requires detailed comparison with cosmological 
simulations; this is beyond the scope of the current paper. 

Figure 2 also shows the detection of an intervening Mgll absorption system. The 
Mgll doublet A2796.4 + 2803.5 is detected at wavelengths 9166.7 A and 9190.3 A in the 
high-resolution spectrum; the redshift of this system is z a b s = 2.278. The rest frame 
equivalent widths of the doublet lines are 2.43 and 1.90 A, respectively. This system is 
very similar to the one detected in SDSSp J033829. 31+002156. 3 (z a b s = 2.304, |Songaila~et 



al. 1999 ). It is possible that SDSS 1044-0125 is amplified by lensing from this intervening 
absorber. However, we saw in § 2 that this object is unresolved under 0.4" seeing in the K 
band. 



4.2. Number Density of Very High Redshift Quasars 

The total area of SDSS imaging data that we have searched for high-redshift quasars 
thus far is of order 600 deg 2 . All that satisfy z* < 19.3 and i* — z* > 2.2 in this 600 
deg 2 region have been observed spectroscopically. Only SDSS 1044-0125 is identified as 
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a high-redshift quasar; the remaining objects are L and T dwarfs. Using the luminosity 
function and redshift dependence of |Schmidt, Schneider fc Gunn (1995)1 (f° r the Q = 1 
model), extrapolating it to higher redshifts and assuming /„ oc z/ -0 - 5 , we predict that in a 
total area of 600 deg 2 , for z > 5.65 (i'-dropout objects ), there should be 1.5 quasars with 
^1450 < —27.0 and 1.1 quasars with M U50 < —27.2. For z > 5.8, the extrapolation predicts 
1.1 and 0.8 quasars for M 1450 < —27.0 and —27.2, respectively. This assumes that our 
selection efficiency is 100%. The Schmidt-Schneider-Gunn luminosity function is derived 
using objects with 2.7 < z < 4.7 and —27.5 < Mb < —25.5. Although it is difficult to draw 
any reliable conclusion from the observation of a single high-redshift quasar, the discovery 
of SDSS 1044-0025 is consistent with the expectations from this rather large extrapolation 
from lower redshift results. Assuming that the same luminosity function holds at even 
higher redshift, the SDSS will be able to discover one quasar at z ^ 6, z* <^ 19 in every 
1500 deg 2 of the survey. 

SDSS 1044-0125 is a very luminous quasar. Assuming that (1) its bolometric luminosity 
equals the Eddington luminosity, Lboi = -^Edd = 1-5 x 1O 38 (Mbh/M ) erg s _1 ; (2) its 
intrinsic continuum spectrum is the same as the mean spectral template of Elvis ct| 



al. (1994")| , and (3) neither beaming nor lensing affects the observed flux, we find a black 



hole mass of M BH = 3.4 x 1O 9 M in the A-model, or M BH = 2.7 x 1O 9 M in the ft = 1 
model. In either case, the implied black hole mass is quite large, similar to that of the black 
hole at the center of the nearby giant elliptical galaxy M87 (Harms et al. 1994, Macchetto 
et al. 1997). If the quasar were radiating below the Eddington limit, the inferred black hole 
mass would be even higher. In the Elvis et al. (1994]| template, ~ 1% of the bolometric 



luminosity is emitted in the observed z' band. If this fraction is larger for SDSS 1044-0125, 
the implied black hole mass would be reduced. Note that the universe was less than 1 Gyr 
old at this redshift, while the Eddington time scale, the e-folding time for the growing 
of a black hole shining at the Eddington luminosity, is 4 x 10 7 (e/0.1) yr, where e is the 
radiative efficiency for the accretion. If the black hole started accreting with an initial 
mass of ~ 1O 3 M with 10% efficiency, the seed black hole would have to form and begin 
accreting at redshift well beyond 10 in order to grow to 3 x 10 9 M Q at z = 5.80 (see also 
Turner 199 1| ) . Forming such a massive black hole in such a short time is a remarkable feat. 



The observations of high-redshift quasars can be used to constrain the formation epoch of 
the first star clusters and the fueling process of early black holes. 

How likely is it to find a quasar like SDSS 1044-0125 in popular cold dark matter 
cosmological models? We use the following simple model to estimate the abundance of 
high-redshift quasars (see also |Haiman fc Hui 2000| ) : [Magorrian et al. (1998)| have found 



a correlation between central black hole mass and bulge mass, MeH/^buige = 6 x 10 3 in 
nearby galaxies. If this correlation holds at high redshift, this would imply a bulge mass of 
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5.7 x 10 11 M (A-model) or 4.5 x 10 11 M (Q = 1) for the host galaxy of SDSS 1044-0125, 
and a lower limit of 5.7 x 10 12 M (A-model) or 4.5 x 10 12 M Q (ft = 1) for its dark halo, 
assuming M ha i /M bu i ge > Q D M/^b ~ 10. 

The comoving abundance of dark matter halos at this epoch is very sensitive to this 
halo mass, and can be estimated by means of the Press fc Schechter (1974)1 formalism. In a 



A-CDM model, assuming <jg = 0.87 and an untilted primordial power spectrum, for parent 
halo masses of 10 12 M , 6 x 10 12 M , and 10 13 M , we expect 50,000, 60, and 4 candidate 
halos respectively within the survey volume in a redshift window Az = 1. The duty cycle 
of quasar activity is poorly known, but is certainly much less than unity ( [Haiman fc Hui| 
2000|). Given these uncertainties, this model is not contradictory to our discovery of a single 



quasar like SDSS 1044-0125 in the 600 deg survey area. 

Note that this model does not address the physical process by which the massive black 
hole formed. The assumptions we made about MeH/^buige and the lifetime of the quasars 
are completely untested at high redshift. Indeed, Rix et al. (1999) argue that the host 
galaxies of z ~ 2 quasars are appreciably less luminous than the universal M B H/^bui g e 
hypothesis would imply. Because of its high luminosity, SDSS 1044-0125 likely probes the 
exponential, high-mass tail of the underlying dark halo distribution, making predictions 
of the expected number counts also sensitive to cosmological parameters, especially the 
normalization of the power spectrum (i.e. <r 8 ). In principle, the SDSS survey will be able 
to probe quasars ~ 1 mag fainter than SDSS 1044-0125, and could reveal tens of additional 
sources at z ~ 6. The detection of these objects will yield strong constraints on cosmological 
models for the formation and evolution of quasars at very high redshifts. 
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Table 1. SDSS and K-band Photometry of SDSSp J10:44:33.04 -01:25:02.2 



u* 


9* 


r* 


i* 


z* 


K 


22.91 


23.91 


25.13 


21.81 


19.23 


17.02 


± 0.53 


± 0.49 


± 0.62 


± 0.19 


± 0.07 


± 0.04 



The SDSS photometry (u*,g*,r*,i*,z*) is reported in terms of asinh magnitudes on the 
AB system. The asinh magnitude system is defined in Lupton, Gunn & Szalay (1999); it 
becomes a linear scale in flux when the absolute value of the signal-to-noise ratio is less than 
about 5. In this system, zero flux corresponds to 24.24, 24.91, 24.53, 23.89, and 22.47, in u*, 
g*, r*, i*, and z*, respectively; larger magnitudes refer to negative flux values. The K band 
photometry is on the Vega-based system. E(B — V) in this direction is 0.054, from Schlegel 
et al. 1998 



Table 2. Emission Line Properties 



Line 


A (A) 


redshift 


EW (A, rest frame) 


OI+SiII 1302 


8858.5 ± 3.8 


5.802 ±0.003 


3.4 ± 0.5 


SilV+OIV] 1400 


9507.3 ± 3.7 


5.789 ±0.003 


7.0 ± 0.4 
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Figure 1. Finding chart for SDSS 1044-0125 (discovery image from the SDSS). The field is 
160" on a side. The field is given in both the %' and z' bands (54.1s exposure time). 
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Figure 2. Optical spectrum of SDSS 1044-0125 observed with Keckll/ESI. The total 
exposure time is 3600s. The spectrum is smoothed to 4 A/pixel. The spectral resolution is 
~ 8 A at A = 9000A. 
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Figure 3. ESI Echellete spectrum over the range 7300 - 8800 A. The spectrum has been 
smoothed to a resolution of 2 A. The dotted line shows a f u oc z/~°' 5 continuum, normalized 
to the region 8650 - 8730 A. The dashed line indicates the region with minimum absorption 
(r ~ 0.4, see §3.1). The lower panel shows the noise level. 



